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Phosphorylation as a Tool To Modulate Aggregation 
Propensity and To Predict Fibril Architecture 
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Despite the importance of post-translational modifications in 
controlling the solubility and conformational properties of pro- 
teins and peptides, precisely how the aggregation propensity 
of different peptide sequences is modulated by chemical 
modification remains unclear. Here we have investigated the 
effect of phosphorylation on the aggregation propensity of a 
13-residue synthetic peptide incorporating one or more phos- 
phate groups at seven different sites at various pH values. 
Fibril formation was shown to be inhibited when a single 
phosphate group was introduced at all seven locations in the 
peptide sequence at pH 7.5, when the phosphate group is 
fully charged. By contrast, when the same peptides were ana- 
lysed at pH 1.1, when the phosphate is fully protonated, fibrils 
from all seven peptide sequences form rapidly. At intermediate 
pH values (pH 3.6) when the phosphate group is mono-anion- 
ic, the aggregation propensity of the peptides was found to 



Introduction 

The aggregation of proteins and peptides into amyloid fibrils is 
associated with a number of life-threatening human dis- 
orders^ 11 including type II diabetes, [2] and Alzheimer's [3] and Par- 
kinson's diseases. [4] The monomeric precursors associated with 
these diseases range from small peptides (A|3, amylin) to na- 
tively unfolded (a-synuclein) and natively folded proteins (|3 2 - 
microglobulin ((3 2 M), lysozyme). Despite their lack of sequence 
identity or structural similarity, the fibrils generated from differ- 
ent proteins share a common cross-|3 structural motif charac- 
terised by hydrogen-bonded (3-strands arranged perpendicular 
to the fibril axis to form |3-sheets, which can subsequently 
stack together to generate a regular assembly. [5] 

Atomic details of the architecture of amyloid fibrils have 
been obtained by using a range of techniques, including solid- 
state NMR spectroscopy (ssNMR), [6] cryo-electron microscopy 
(cryo-EM), [7] electron paramagnetic resonance (EPR), [8] X-ray 
fibre diffraction/ 91 hydrogen/deuterium (H/D) exchange/ 101 and 
proline-scanning mutagenesis/ 111 Despite these developments, 
it is still a challenge to obtain atomistic details about amyloid 
structures. To meet these challenges, recent studies have fo- 
cused on the use of small peptides as model systems to gain 
atomic details on the organisation of the |3-strands and |3- 
sheets within the fibril. X-ray crystallographic analysis of fibrils 
formed from a range of short amyloidogenic peptide frag- 
ments has highlighted a limiting set of eight possible arrange- 
ments of (3-strands and (3-sheets within the cross-p structure of 
amyloid, referred to as steric zippers because of the tight inter- 



be highly dependent on the position of the phosphate group 
in the peptide sequence. Using this information, combined 
with molecular dynamics (MD) simulations of the peptide 
sequence, we provide evidence consistent with the peptide 
forming amyloid fibrils with a class 7 architecture. The results 
highlight the potential utility of phosphorylation as a method 
of reversibly controlling the aggregation kinetics of peptide 
sequences both during and after synthesis. Moreover, by ex- 
ploiting the ability of the phosphate group to adopt different 
charge states as a function of pH, and combining experimental 
insights with atomistic information calculated from MD simula- 
tions as pH is varied, we show how the resulting information 
can be used to predict fibril structures consistent with both 
datasets, and use these to rationalise their sensitivity of fibrilla- 
tion kinetics both to the location of the phosphate group and 
its charge state. 



meshing of the p-sheet side chains to form a dry interface 
within the protofilament structure/ 121 Although this approach 
has proved successful in determining the structures of a large 
number of proteins or peptides, they rely on the ability to crys- 
tallise relevant small fragments. Hence, a simple, rapid and re- 
versible chemical approach to gain insight into fibril structure 
and stability is warranted. 

Phosphorylation has been used to modulate the aggrega- 
tion propensity of proteins and peptides, including tau, [13] a- 
synuclein/ 141 and peptide model systems/ 151 The amyloidogenic 
potential of a-synuclein, for example, can be suppressed, in 
vitro, at neutral pH by inserting a phosphate group at Ser129; 
however, mutating the serine residue into an aspartate or a 
glutamate does not mimic the effect of the phosphate group, 
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despite their structural and electrostatic similarities^ 141 Phos- 
phorylation has also been used to modulate the aggregation 
propensity of a coiled coil peptide model by enzymatically trig- 
gering a structural switch leading to amyloid formation. [15bc] 
However, these studies have provided only a phenomenologi- 
cal description of the effect of the addition of a phosphate 
moiety, and they do not consider in detail either pH or posi- 
tional dependence or an atomistic description of the observed 
effects. 

In contrast, in this study, we report for the first time, the se- 
quence- and pH-dependence of the effect of phosphorylation 
on aggregation using a small peptide of the human protein 
|3 2 M corresponding to residues 59 to 71, and by combining ex- 
periment and simulation we provide atomistic models for the 
origin of the effects observed. Although there is no evidence 
that the full-length protein is phosphorylated, in vivo, this 13- 
residue sequence (DWSFYLLYYTEFT, known as strand E) [16] 
forms an ideal model system for the studies presented, as six 
out of the 13 residues can be phosphorylated reversibly, in 
vitro, and the peptide has been shown to be highly aggre- 
gation-prone when using theoretical methods [17] and in vitro 
aggregation studies. [16,18] The effect of inserting a phosphate 
moiety at a specific position within the sequence on fibril for- 
mation was then investigated by using thioflavin T (ThT) fluo- 
rescence and negative stain transmission electron microscopy 
(TEM) at pH values corresponding to the different protonation 
states of the phosphate group. Protofilament structures of the 
fibrils formed by strand E and its phosphorylated variants con- 
sistent with the experimental data were then constructed and 
their stabilities as a function of pH were monitored by using 
MD simulations. By combining the experimental and simula- 
tion data, we derive a model for the fibrillar architecture of this 
peptide able to satisfy all the data obtained that contains an 
antiparallel, shifted organisation of the |3-strands organised 
within a fibril of class 7 [12] architecture. 

Results and Discussion 

Synthesis of peptides 

In order to investigate the effect of phosphorylation on the ag- 
gregation propensity of the amyloidogenic sequence of |3 2 M 
corresponding to residue 59 to 71 (strand E), a series of tride- 
capeptides containing zero, one or two phosphoamino acid 
residues was synthesised by using Fmoc solid-phase synthe- 
sis. [19] The C terminus was prepared as the amide and the N ter- 
minus was capped with an acetyl group. Significant optimisa- 
tion of the synthesis was required as high loading resins led to 
incomplete amino acid coupling when the core amyloidogenic 
sequence "FYLLYY" was reached. Synthesis of peptides with a 
high tendency to aggregate and form intermolecular |3-sheet- 
like structures often leads to incomplete amino acid coupling, 
and therefore, poor overall yields. [19b,20] Strand E has previously 
been demonstrated to be highly aggregation-prone, [16] which 
makes it a particularly "difficult sequence" to synthesise. 

Low yields and aggregation of the peptide sequence within 
the resin matrix were prevented by employing a low loading 
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resin and inserting the |3-sheet breaker 2-hydroxy-4-methoxy- 
benzyl (Hmb) protecting group prior to the amyloidogenic 
sequence "FYLLYY" at position 12. [21] Use of a highly acid labile 
resin (NovaSyn TG Sieber resin) also allowed the fully depro- 
tected peptide variants to be obtained by using a two-step 
cleavage procedure (Figure S1 in the Supporting Information). 
The peptides were first released from the solid support in their 
fully protected forms with a low concentration (2% v/v) of tri- 
fluoroacetic acid (TFA). During this stage the aggregation-dis- 
rupting effect of the Hmb protecting group is retained, and 
prevents aggregation of the peptide variants within the matrix 
of the solid support. After release of the peptide from the sup- 
port, global deprotection was achieved with a high concentra- 
tion (94% v/v) of TFA and appropriate cation scavengers. The 
series of peptides shown in Table 1 were finally purified by 
HPLC and shown to be the correct mass by ESI-MS. 



Table 1. 


Peptides synthesised in this investigation. 131 
Peptide target 


Yield [%] [b] 


control 


DWSFYLLYYTEFT 


23 


p13T 


DWSFYLLYYTEFpT 


10 


plOT 


DWSFYLLYYpTEFT 


8 


p9Y 


DWSFYLLYpYTEFT 


18 


p8Y 


DWSFYLLpYYTEFT 


8 


P 5Y 


DWSFpYLLYYTEFT 


8 


P3S 


DWpSFYLLYYTEFT 


12 


p3Sp10T 


DWpSFYLLYYpTEFT 


9 


[a] Phosphopeptide variants are labelled pX, where X denotes the phos- 
phorylated amino acid side chain. All peptides were N-terminally acetylat- 
ed and C-terminally amidated. [b] Quoted yields are of pure peptide ob- 
tained after purification by preparative HPLC. 



Aggregation of nonphosphorylated strand E (control) 

The aggregation kinetics of control were measured at pH 1.1, 
3.6 and 7.5 by using ThT fluorescence and TEM as a probe of 
fibrillogenesis (Figure 1). These pH values were chosen to 
study the effect of three charge states of the phosphoamino 
acid moiety on the aggregation propensity of the peptides, 
ranging from 0, at acidic pH, to —2, at neutral pH. Aggregation 
of control was initiated by addition of a buffer solution con- 
taining ThT to a peptide stock solution in dimethylsulfoxide 
(DMSO) and the ThT fluorescence was measured in real-time in 
a 96-well plate format (Figure 1 A). The presence of fibrils was 
subsequently confirmed by TEM (Figure 1 B). 

Analysis of the aggregation kinetics of control showed that 
the peptide formed fibrils rapidly at all pH values studied. The 
rate of fibril formation was then determined by using 40 repli- 
cate curves by extracting an average half-life value (f 50 ) with 
the process described in Figure 2. A f 50 value was extracted for 
each of the 40 normalised curves and these values were plot- 
ted on a histogram, which was then fitted to a Gaussian distri- 
bution to generate an average f 50 value for each peptide under 
each set of conditions (Figure S2 and Table S1 in the Support- 
ing Information). 
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Figure 1. Aggregation of the non-phosphorylated peptide control at pH 1.1, 
3.6 and 7.5 at 30 °C. A) Spontaneous fibril growth monitored by ThT fluores- 
cence, and B) negative-stain TEM images of control after incubation for 
4 days at the different pH values (scale bar: 100 nm). The kinetic curves 
shown in triplicate are representative of 40 replicates for each peptide ob- 
tained at each pH value. The ThT curves are normalised to the ThT fluores- 
cence of samples at pH 1.1. 
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The native sequence of strand E contains two negatively 
chargeable residues, aspartate and glutamate, at positions 1 
and 11, respectively. Assuming pK a values of approximately 
four for these groups, the net charge of control varies from 
0 at pH 1.1 to —2 at pH 7.5, which correlates well with the 
increase of the f 50 value between these pH values 
(0.367(±0.003) s at pH 1.1, 0.484(±0.006) s at 3.6, and 
10.59(±0.18) s at pH 7.5). The presence of the two negative 
charges at the aspartate and glutamate residues thus retard 
significantly the rate of aggregation of the peptide, but are in- 
sufficient to inhibit fibril formation completely, even when fully 
charged at neutral pH (Figure 1 A and B). 

The effect of the phosphate group is electrostatic rather 
than steric 

To assess the effect of a single phosphate moiety on the ag- 
gregation kinetics of strand E, the phosphopeptide variant p9Y 
containing a single phosphotyrosine at position 9 was pre- 
pared and its aggregation kinetics were measured (Figure 3). 

Whereas control formed fibrils rapidly at all pH values stud- 
ied, the insertion of a phosphate group at position 9 into the 
sequence of strand E had a dramatic effect on assembly. While 
p9Y aggregated rapidly at pH 1.1 (f 50 = 2.90(±0.02) s) and 
formed long straight fibrils as confirmed by TEM (Figure 3 B), 
the aggregation kinetics were retarded at pH 3.6 (f 50 = 

74.77(±0.77) s) and assembly 
was completely inhibited when 
the pH was increased to 7.5 (Fig- 
ure 3 A and B). The ability of p9Y 
to form fibrils rapidly at pH 1.1 
indicates that incorporation of a 
phosphate moiety at this site 
can be accommodated within 
the fibrillar architecture. Rather 
than being sterically determined, 
the inhibition of fibrillar assem- 
bly only occurs when the phos- 
photyrosine is in a negatively 
charged state. 
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Figure 2. Process used to extract an average f 50 value for control at pH 1.1. A) Raw data showing the 40 replicates 
after subtracting the DMSO blank, B) normalised data showing the 40 replicates after normalizing each curve 
against the average value of its plateau, C) extraction of a f 50 value from a single normalised curve, and D) histo- 
gram of f 50 values fitted with a Gaussian (black line) for all 40 replicate curves obtained. An average f 50 value was 
calculated from the fit to the Gaussian function (Table SI in the Supporting Information). 



The position of the phosphate 
group is an important modula- 
tor of aggregation propensity 

With these data in hand, we 
next investigated whether phos- 
phorylation can be used as a 
general strategy for modulation 
of aggregation and whether the 
effect is sequence dependent. 
Accordingly, the phosphate 
group was introduced at seven 
different positions within the 
amino acid sequence of strand E 
(Table 1) and the aggregation ki- 
netics of the different peptide 
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Figure 3. Fibril formation of the phosphopeptide variant p9Yat pH 1.1, 3.6 
and 7.5 at 30 °C. A) Spontaneous fibril growth monitored by ThT fluores- 
cence, and B) negative stain TEM images after incubation for 4 days (scale 
bar: 100 nm). The kinetic curves shown in triplicate are representative of 40 
replicates obtained at each pH value. 



variants were measured at pH 1.1, 3.6 and 7.5, as described for 
control (Figures 4 and 5). At pH 1.1, all the phosphopeptides 
aggregate rapidly to form amyloid-like fibrils, as measured by 
ThT fluorescence (Figure 4A) and confirmed by negative stain 
TEM (Figure S3 in the Supporting Information). For all peptides 
containing a single phosphate group, the aggregation kinetics 
at pH 1.1 have f 50 values below 10 s. For the variant containing 
two phosphate groups, p3S10T, the kinetic curves show a 
wide range of aggregation rates, ruling out determination of a 
single f 50 value. Similar to the behaviour of p9Y, the aggrega- 
tion kinetics of all other phosphopeptides was completely sup- 
pressed at pH 7.5 (Figure 4C). 

At pH 3.6, aggregation was highly dependent on the se- 
quence, with each phosphopeptide variant exhibiting a differ- 
ent kinetic behaviour (Figure 4 B). Whereas four of the mono- 
phosphorylated peptides (p13T, plOT, p9Y and p3S) form fi- 
brils at pH 3.6, two variants (p5Y and p8Y) were unable to 
form fibrils at this pH even after incubation for 4 days. Peptides 
p3S and pi 0T aggregated at pH 3.6 (Figure 5); however, the 
equivalent double phosphopeptide p3S10T did not form any 
fibrils under these conditions (data not shown). 

Inhibition of amyloid formation provides constraints on the 
fibril architecture 

The sequence-dependent effect observed at pH 3.6 for the 
phosphopeptide variants described in Table 1 led us to consid- 
er whether the aggregation rates of the different peptides 
could provide information about the possible fibrillar architec- 
ture adopted by strand E. Assuming that electrostatic repulsion 
between adjacent phosphate moieties modulates fibril assem- 
bly, and that the different peptide sequences assemble into a 



A) 



3.2 
[ 2.8 

£ 2.0 
| 1.6 
I- 1.2 

f 0-8 
° 0.4 
0.0 




012345678 

f/min — ► 




B) 



| 2.4* 



2.4 
2.0 
1.6 
1.2 
0.8 



~ 0.4 




0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 

f/min — ► 




Figure 4. ThT kinetic curves of seven peptide variants of strand E at: 

A) pH 1.1, B) pH 3.6, and C) pH 7.5. Each peptide is coded by colour: control, 

, plOT, p9Y, p8Y, p5Y, p3S and p3S10T. Representative negative-stain 
TEM images are shown after incubation of control, , and p5Y for 4 days 
at 30 °C (scale bar: 100 nm). Electron micrographs for each peptide at each 
pH value are shown in Figure S3 in the Supporting Information. The kinetic 
curves shown in triplicate are representative of 40 replicates for each pep- 
tide obtained at each pH value. An inset expansion is shown for pH 7.5, 
showing the first 6 min of assembly. Note the time axis in graphs (A), (B) 
and (C) are different. 
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Figure 5. Histogram of average f 50 values of aggregation of strand E pep- 
tides containing a single phosphorylated amino acid at pH 1.1, 3.6 and 7.5. 
Peptides showing no visible fibrils at the end of the incubation time (4 days) 
are shown with a f 50 value of > 4 days. The f 50 values are plotted on a loga- 
rithmic scale. Standard deviations depicted as error bars are too small to be 
seen on this scale (Table SI in the Supporting Information). 



common fibril morphology, as is frequently (but not always) 
found for close sequence homologues, [22] the different assem- 
bly properties of the peptides at pH 3.6 (under which condi- 
tions their assembly rates are most different) can be exploited 
to derive models of the fibril architectures most likely to be 
formed. 

Eisenberg et al. have defined a limiting set of eight possible 
arrangements of peptides in a cross-|3 structure. These struc- 
tures are referred to as "steric zippers" because the side chains 
of residues in adjacent |3-sheets intermesh tightly like the teeth 
of a zipper to form a dry interfaced 121 These structures are for- 
mally defined by the symmetry relationships in the organisa- 
tion of the p-strands and |3-sheets with respect to one another. 
We first considered the first level of the steric zipper architec- 
ture — the arrangement of p-strands within a p-sheet — to de- 
termine whether the fibrillation kinetic data for different phos- 
phovariants of strand E could support, or eliminate, certain 
classes. The most simplistic representation of a steric zipper 
consists of a pair of p-strands that can be organised in one of 
three arrangements. While only one parallel arrangement is 
possible, a p-strand can generate two antiparallel sheet struc- 
tures, depending on whether the p-strands are related by two- 
fold rotational (C 2 ) symmetry around the backbone axis or per- 
pendicular to the backbone axis. The alternative operations 
modulate the orientation of side chains that project from the 
p-strands and also lead to a shifted or in-register organisation 
of the p-strands (Figure 6). [23] 

Three simplified strand-strand schematic models of control 
were constructed in these three strand architectures with the 
tyrosine residues coloured according to the position of the 
phosphate groups in p5Y, p8Y and p9Y (Figure 6). These pep- 
tide variants were chosen because they show very different 
rates of aggregation at pH 3.6, with both p5Y and p8Y (orange 




Figure 6. Schematic representation of the possible strand-strand interac- 
tions within a (3-sheet for strand E with the p-strands: A) parallel, in-register, 
B) antiparallel, in-register, and C) antiparallel, shifted arrangement. The col- 
oured amino acids represent the tyrosine residue in which the phosphate 
groups would be present: p9Y (green), p8Y (orange) and p5Y (purple). The 
red arrow shows the long axis of the fibrillar array. 

and purple) being incapable of forming fibrils at this pH, while 
p9Y (green) aggregates rapidly. 

Assuming that modulation of the aggregation kinetics of the 
different phosphopeptides results from electrostatic repulsion 
between adjacent phosphate groups, the proximity of the 
phosphotyrosine residues was used to determine whether 
peptide E is most likely to adopt parallel or antiparallel p- 
strands in the fibril structure. Peptides p5Y and p8Y do not ag- 
gregate at pH 3.6, whereas p9Y is able to form fibrils at this 
pH. Using these criteria and experimental observations, a paral- 
lel, in-register organisation of p-strands (Figure 6A) can be 
ruled out since in this arrangement of p-strands each phos- 
phorylated amino acid is positioned equivalently within the p- 
sheet. Similarly, an antiparallel, in-register arrangement of 
strands within the p-sheet (Figure 6 B) cannot account for the 
different aggregation properties of peptides p5Y and p9Y as 
the tyrosine residues at positions 5 and 9 are equally spaced. 
The only p-sheet arrangement that can account for the differ- 
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ent behaviours experimentally observed for p5Y, p8Y and p9Y, 
therefore, is the antiparallel shifted organisation within a 13- 
sheet (Figure 6C). This model accounts fully for the similar ag- 
gregation kinetics of p5Y and p8Y since the phosphotyrosine 
groups are spaced identically in this p-sheet structure. By con- 
trast, the separation of the tyrosine residues bearing the phos- 
phate group in p9Y is greater in this fibril structure accounting 
for the ability of p9Y to form fibrils at pH 3.6. This structure 
also allows the maximum jt-jt stacking interactions in the core 
sequence, and the alignment of the hydrophobic leucine resi- 
dues — features presumed to play an important role in deter- 
mining fibril stability. [24] 

Of the eight classes of steric zippers described by Eisenberg 
et al., only two classes possess the p-strand symmetry of the 
shifted organisation described in Figure 6C: the class 7 or 8 ar- 
chitectures.^ We, therefore, considered if the data presented 
above can allow discrimination between the different (3-sheet 
packing modes within the fibril architectures described by 
these alternative classes. While the arrangement of a pair of (3- 
strands can allow certain models to be ruled in or out as dis- 
cussed above, these simplistic strand-strand schematics do 
not account for the contacts occurring at the p-sheet interface. 
Models of aggregates equivalent to a pair of interdigitating p- 
sheets of eight p-strands were therefore built in silico for fibrils 
with class 7 and 8 architectures with the nonphosphorylated 
peptide (control; Figure 7), assuming an antiparallel shifted p- 
strand organisation within each p-sheet. The difference be- 
tween these two steric zippers relates to the relative orienta- 
tion of the faces of the two identical p-sheets. In order to de- 




Figure 7. The 8x2 steric zippers of strand E control showing the intersheet 
packing for: A) class 7, and B) class 8, as described by Eisenberg et al. [12] The 
coloured amino acids represent the tyrosine residues that are phosphorylat- 
ed in p9Y (green), p8Y (orange) and p5Y (purple). The long axis of the fibril- 
lar array is shown as a red dot, the fibrillar axis goes into the plane of the 
paper. 
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termine which class of steric zippers is the more likely arrange- 
ment for strand E fibrils, an average distance between the Cp 
atoms of the tyrosine residues at positions 5, 8 and 9 was cal- 
culated (Table 2). 



Table 2. Average C p -C p distance between tyrosine residues at positions 5, 


8 and 9 in a steric zipper of class 7 and 8. 




Peptide Class 7 C p -C p [A] 


Class 8 C p -C p [A] 


5Y 13.7 


11.6 


8Y 13.6 


10.4 


9Y 20.5 


10.6 



The experimental results obtained at pH 3.6 show that pep- 
tides p5Y and p8Y are unable to form fibrils (Figure 5), while 
p9Y aggregates rapidly at this pH; this suggests that the phos- 
photyrosine groups of p5Y and p8Y must be in closer proximi- 
ty than that of p9Y in the intersheet packing in the fibril struc- 
ture. In a class 8 architecture (Figure 7 B), the average Cp-Cp 
distance between the tyrosine residues of 5Y, 8Yand 9Y (corre- 
sponding to the position of the three tyrosine residues where 
the phosphate groups have been introduced) is comparable 
(Table 2). By contrast, in a class 7 architecture, 5Y and 8Y have 
comparable Cp-Cp distances, while the tyrosine residues of p9Y 
are significantly further apart. The class 7 architecture is there- 
fore more consistent with the aggregation propensities of the 
peptides phosphorylated at these residues. Analysis of these 
two models thus suggests that a steric zipper with a class 7 
architecture is the only organisation consistent with the experi- 
mental results obtained (Figure 7A). [12] 

A class 7 architecture as a viable polymorph of strand E 

If a polymorph corresponding to a class 7 architecture is pref- 
erentially formed across the peptides and pH range studied, 
the stability of the models built for the phosphopeptide var- 
iants p5Y, p8Y and p9Y with a class 7 architecture should cor- 
relate directly with the aggregation behaviours observed ex- 
perimentally. To determine whether this is the case and to pro- 
vide a deeper insight into the inter-atomic interactions within 
fibrils containing phosphotyrosine moieties at different sites 
within a common class 7 architecture, in silico models of fibrils 
containing eight p-strands and two p-sheets with charged and 
protonated phosphotyrosine residue at positions 5, 8 or 9 
were built. Using atomistic MD simulation, we then monitored 
whether the fibril remained as a stacked pair of ordered p- 
sheets as would be expected for stable architectures. The sim- 
ulations are used as a design tool based on the assertion that 
if the electrostatic repulsion between charged phosphate moi- 
eties within the fibril is sufficiently large to prevent peptide ag- 
gregation experimentally, the originally ordered fibrils built in 
silico will collapse into a disordered state during the MD calcu- 
lations. MD trajectories 10 ns in length in explicit solvent were 
therefore obtained for the structures bearing either fully pro- 
tonated, mono- or di-anionic phosphotyrosine groups at each 
site. Figure 8 shows snapshots of the structures at 0, 5 and 
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10 ns of the MD simulations in 
which the solvent distribution 
around the leucine and the 
phosphotyrosine residues of 
p5Y, p8Y and p9Y are highlight- 
ed. The results are striking, and 
suggest that the solvent and 
counter-ion environment plays 
an important role in determining 
the behaviour of the fibrils deco- 
rated with phosphotyrosine 
groups carrying differing 
amounts of charge, as discussed 
in detail below. 

With a fully protonated phos- 
photyrosine group, which is 
electrostatically neutral (Fig- 
ure 8A), the steric zippers of 
p5Y, p8Y and p9Y showed little 
or no variation in their architec- 
tures throughout the simulation, 
and the dry, intersheet interface 
remained solvent-free over the 
duration of the MD. As no water 
molecules penetrated the inter- 
face the side chains of both (3- 
sheets remained dry and tightly 
packed. These results suggest 
that the phosphopeptide var- 
iants p5Y, p8Y and p9Y bearing 
a fully protonated phosphotyro- 
sine group can adopt a class 7 
architecture and accommodate a 
phosphate group inside the 
tight dry interface while remain- 
ing in an ordered fibrillar state. 
Transient hydrogen bonding in- 
teractions between phospho 
groups within the interface 
could be observed in the struc- 
ture of p5Y and to a lesser 
extent in the structure of p8Y. In 
addition, the phophotyrosine 
group in p8Y and p9Y zippers 
showed transient intersheet hy- 



Figure 8. Snapshots of the MD simula- 
tions carried out in explicit solvent for 
the phosphopeptide variants p5Y, p8Y 
and p9Y bearing: A) a fully protonated, 
B) a mono-anionic, and C) a di-anionic 
phosphate group. Snapshots taken at 
f = 0, 5 and 10 ns were prepared with 
VMD. [25] The water molecules (red dots) 
are represented as a surface area within 
10 A of residues Leu and the phospho- 
tyrosine, and the yellow spheres repre- 
sent the sodium ions used to neutralise 
the overall structure. 
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drogen bonding to the hydroxyl group of Y9 and Y8, respec- 
tively, during the simulations. The MD simulations are thus 
consistent with the aggregation behaviour observed experi- 
mentally for p5Y, p8Y and p9Y at pH 1 .1 . 

By replacing the fully protonated phosphotyrosine group 
with its di-anionic form (Figure 8 C), thereby introducing a -2 
electrostatic charge per phosphotyrosine group, the stability of 
the steric zippers of p5Y, p8Y and p9Y was greatly affected. 
During the simulation, water molecules and sodium counter- 
ions penetrated the dry interface, leading to the dissociation 
of each steric zipper. The phosphopeptide variants p5Y, p8Y 
and p9Y bearing a di-anionic phosphate group are therefore 
energetically unstable, such that electrostatic repulsion and/or 
preferential solvation of the negatively charged phosphotyro- 
sine groups results in disassembly of these structures at neu- 
tral pH. These results correlate with the inhibition of the aggre- 
gation of p5Y, p8Y and p9Y observed by ThT fluorescence at 
pH 7.5. 

Similar to the results obtained experimentally, when the 
phosphotyrosine group was replaced with its mono-anionic 
form a sequence-dependent effect was also observed by MD 
simulations (Figure 8 B). While the steric zippers of p5Y and 
p8Y rapidly disassemble as water molecules and sodium coun- 
ter-ions penetrate, the dry interface formed by the two stacked 
p-sheets, the steric zipper of p9Y remained free of water mole- 
cules and sodium counter-ions throughout the simulation. As 
previously shown in Table 2, the phosphotyrosine groups of 
p9Y are much further apart in the intersheet packing than 
those of p5Y and p8Y, with the 
consequent reduction in electro- 
static repulsion within the steric 
zipper of p9Y in this fibril struc- 
ture. Considering that the inter- 
face formed by the tight packing 
of the side chains must be free 
of solvent or counter-ions for a 
steric zipper to be considered 
stable, the structural integrity of 
the fibrillar models built for p5Y, 
p8Y and p9Y can be quantita- 
tively assessed by calculating the 
number of water molecules pen- 
etrating the dry interface over 
the duration of the trajectories 
(Figure 9). 

These data show that the 
models generated for p5Y, p8Y 
and p9Y bearing protonated 
phosphotyrosine groups re- 
mained dry with less than five 
molecules of water penetrating 
within the |3-sheet interface over 
the duration of the MD simula- 
tions. When the phosphotyro- 
sine groups were replaced with 
their mono-anionic equivalent, 
the water penetration appeared 



to be directly related to the position of the phosphotyrosine 
group in the sequence of strand E. While less than five water 
molecules penetrate the |3-sheet interface of p9Y, the struc- 
tures of p5Y and p8Y are rapidly dissolved following the diffu- 
sion of a large number of water molecules (ca. 60 and 20 
water molecules, respectively) within the interface (Figure 9). 
The dissolution effect observed in the MD simulations with a 
mono-anionic phosphotyrosine group is therefore sequence 
dependent and correlates precisely with the aggregation be- 
haviour observed experimentally at pH 3.6 (Figures 4 B and 5). 
With a di-anionic phosphate group, the simulations show a 
rapid solubilisation of the fibrillar models of p5Y, p8Y and p9Y 
(ca. 55, 35 and 35 waters molecules penetrating the dry inter- 
face, respectively) independent of the position of the phos- 
phate group in the sequence of strand E (Figure 9). Taken 
together, therefore, the MD simulations and experimental data 
of variants p5Y, p8Y and p9Y obtained at pH 1.1, 3.6 and 7.5 
imply that strand E adopts a class 7 fibrillar architecture that is 
uniquely able to accommodate a singly charged phosphotyro- 
sine within its dry (3-sheet interface. 

Conclusion 

The ability of post-translational modifications to modulate the 
aggregation propensity of peptides and proteins in vivo is im- 
portant in the progression of a range of amyloid related dis- 
eases. [14,25] In this study we have significantly advanced the un- 
derstanding of how modulation of amyloid formation by post- 




Figure 9. Number of water molecules found within the (3-sheet interface of the class 7 steric zipper of p5Y, p8Y 
and p9Yover the duration of the MD simulations. Snapshots highlighting the position of the water molecules at 
the end of the simulations are shown for some fibrillar models (snapshots of all three peptides with the three dif- 
ferent protonation states are summarised in Figure S4 in the Supporting Information). The phosphopeptide var- 
iants are coloured accordingly to the position of their phosphate groups, p9Y (green), p8Y (orange) and p5Y 
(purple). The different protonation states of the phosphate group are represented as follow: 0 with solid squares, 
-1 with open squares, and -2 with crosses. 
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translational modifications can be understood at the atomistic 
level. By combining experimental and theoretical studies of a 
model peptide we have revealed new features of how phos- 
phorylation can modulate the aggregation of a peptide se- 
quence. Examining the influence of both position and charge 
of the modified residues on the aggregation rates reveals that 
phosphorylated residues can be sterically accommodated 
within the steric zippers of cross-p assemblies, but that as the 
residues become highly charged aggregation is strongly inhib- 
ited. At intermediate charge states the effect of phosphoryla- 
tion becomes sequence dependent. The sequence-specific 
effects have enabled us to provide predictive insights into the 
possible architecture of the fibrils. By combining the experi- 
mental data with atomistic MD simulation of possible poly- 
morphs we gain a consistent picture in which experimental 
aggregation propensity and measure of fibril stability in silico 
suggest a class 7 steric zipper for the fibril architecture of this 
peptide sequence. 

The results presented show that phosphorylation provides a 
powerful and tunable control of aggregation propensity that 
can be exploited as a switch able to turn on/off the aggrega- 
tion of self-assembling materials. Accordingly, assembly can be 
triggered by simple pH switch, or by removal of the phosphate 
group. [14,15] Indeed, dephosphorylation of p9Y was efficiently 
performed in the presence of alkaline phosphatase at pH 8.9 
to generate dephosphorylated peptide E, which aggregated 
with a similar f 50 to the structurally identical, synthetic equiva- 
lent (data not shown). One of the most fascinating features of 
amyloid is the apparent diversity of structures available to a 
single polypeptide sequence within a generic cross-p architec- 
ture [12,26] a bj|j t y t0 identify which sequence can assemble 
into which polymorph is a major challenge that is crucial for 
the development of molecular understanding of amyloid dis- 
ease and to exploit amyloid as a functional nanomaterial. In 
addition to its practical utility, the data presented suggest that 
phosphoamino acid scanning might be a useful method for 
the establishment of putative fibrillar architectures and poly- 
morph distinction, which can then inform the design of meth- 
ods of structural confirmation. Currently, generation of such 
atomistic descriptions of amyloid is a major undertaking, typi- 
cally involving MAS-NMR spectroscopy and cryo-EM. High reso- 
lution structures of amyloid that are consistent with both long 
and short range information from such experiments, however, 
have been achieved for very few systems to date. [27] In this 
study we have shown that simple chemical modification in 
combination with simulation can be used to develop fibril 
models efficiently and rapidly and will be of wide utility in pro- 
viding rapid development of structural models for other as- 
semblies of other peptides. The simulations also suggest 
mechanisms by which the change in charge state of the phos- 
phorylated residue modulates aggregation. Rather than elec- 
trostatic repulsions simply forcing aggregates apart, charge 
density in the interior of the steric zipper promotes the inva- 
sion of ions and water into the intersheet space promoting dis- 
solution of the assembled architecture. 



Experimental Section 

Peptide synthesis: Peptide variants were synthesised manually by 
using standard solid-phase techniques for Fmoc A/a-protected 
amino acids on a NovaSyn TG Sieber resin (0.20 mmolg" 1 , Nova- 
biochem). [19b] Standard Fmoc-protected and phosphorylated amino 
acids (N-a-Fmoc-O-benzyl-L-phospho-serine, -threonine and -tyro- 
sine, Novabiochem) couplings were carried out with excess 2-(6- 
chloro-1H-benzo-triazole-1-yl)-1,1,3,3-tetramethylammonium hexa- 
fluorophosphate (HCTU, Novabiochem) and A/,A/-diisopropylethyl- 
amine (DIPEA, Aldrich) in A/,A/-dimethylformamide (DMF). The |3- 
sheet breaker A/-a-Fmoc-A/-a-(2-Fmoc-oxy-4-methoxybenzyl)-L-Phe 
(Fmoc-(FmocHmb)Phe-OH, Novabiochem) was inserted at posi- 
tion 12 to prevent on-resin aggregation by using HCTU and DIPEA, 
and the subsequent amino acid was coupled, overnight, with 
excess 2-(7-aza-1 H-benzotriazole-1 -yl)-1 ,1 ,3,3-tetramethylammoni- 
um hexafluorophosphate (HATU, Novabiochem) and DIPEA. Depro- 
tection of the Fmoc-protecting group was achieved with a solution 
of piperidine (20%) in DMF. Couplings and deprotections were 
confirmed by using the Kaiser test. 

Acetylation of the N terminus was achieved with excess acetic an- 
hydride and DIPEA in DMF. Cleavage and side-chain deprotection 
was accomplished in two steps. The peptide was first released 
from the resin with a solution of DCM/TFA (98:2) and the resulting 
solid was finally deprotected in solution with a cleavage cocktail of 
TFA/EDT/H 2 0/PhOH (94:2:2:2), precipitated in cold diethyl ether 
and lyophilised from AcOH (0.1 %). 

The crude peptides were purified by preparative reverse phase 
HPLC on a Gilson instrument with a Phenomenex Jupiter Proteo 
column (10|tim, 90 A, 250x21.2 mm, 20ml_min" 1 ) and eluted with 
a linear gradient of 5% to 25% B, where B is acetonitrile + 0.1 % 
ammonia and A is water+10% B. Analytical traces and accurate 
mass data were obtained to confirm purity. 

Measurement of aggregation by using ThT fluorescence: All 

fluorescence experiments were performed at 30 °C on a Perkin- 
Elmer EnVision 2103 Multilabel reader (excitation filter: photomet- 
ric 440, band width: 8nm; emission filter: FITC 485, band width: 
14 nm; general beam splitter in the optical path; ten flashes per 
well), by using a 96-well plate (MicroFluor). Buffers with a final 
ionic strength of 150 mM (A: 25 mM glycine, 47.0 itim NaCI 
(pH 1.1); B) 25 mM sodium phosphate, 25 mM sodium acetate, 
99.7 mM NaCI (pH 3.6); and C) 25 mM sodium phosphate, 83.5 mM 
NaCI (pH 7.5)) were mixed in a 9:1 ratio with ThT (1 mM in water). 
Peptide stock solutions (2 (iL, 1 mM in DMSO) were diluted to a 
final concentration of 20 jlim with ThT buffer (98 |ul). In the nega- 
tive control the peptide was replaced with DMSO. 

At pH 1.1 and 3.6, buffers A and B, respectively, were dispensed di- 
rectly from the plate reader into a 96-well plate, containing either 
the peptide stock solution or DMSO. At pH 7.5, buffer C was dis- 
pensed with a Hamilton Star robotic liquid handler into a 96-well 
plate and the plate was sealed with the optically clear film Plate- 
max (Axygen). A continuous assay was used to monitor fibril for- 
mation (at pH 1.1, 250 measurements at intervals of 2 s; at pH 3.6, 
170 measurements at intervals of 5 s apart for plOT (10 s intervals); 
at pH 7.5, 100 measurements at intervals of 4 min). 

The fluorescence intensity of each sample was corrected for the 
ThT signal of the negative control (Figure S2 in the Supporting In- 
formation). Control experiments with preformed fibrils (incubated 
in neat buffer at a final concentration of 40 jlim for 4 h) demon- 
strated that ThT binding to the preformed fibrils occurs within the 
dead time of less than 0.5 s at all pH values studied. 
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The f 50 determination: The 40 kinetic curves were each normalised 
to the average of their end point (plateau), and because they 
could not be analytically fitted, the time to reach 50% of the maxi- 
mum signal (f 50 ) was numerically calculated for each curve by 
using Equation (1): 

r 50 = r 1+( o.5-y 1 )(£=£) 0) 

where r, represents the time point below f 50 and t 2 the time point 
above; and y 2 are the normalised fluorescence intensities at f, 
and f 2 , respectively. 

For each peptide at each pH, a frequency count was carried out 
for the 40 f 50 values obtained with a total number of bins equal to 
10 and a histogram was plotted. A Gaussian curve was fitted to 
each histogram, and mean f 50 values and standard deviations were 
calculated for each peptide for each set of conditions based on a 
normal distribution. 

Transmission electron microscopy (TEM): Samples were prepared 
in an Eppendorf tube by diluting peptide stock (10|liL; 1 itim in 
DMSO) in buffer (500 jliL). After incubation for 4 days at 30 °C with 
shaking at 200 rpm, the solution (10 ui) was dispensed undiluted 
on to freshly ionised formvar- and carbon-coated TEM grids (Agar) 
and incubated for 1 min. The grids were stained for 30 s with 
uranyl acetate (4% w/v). All images were taken by using a Philips 
CM10 electron microscope operating at 80 keV. 

Model minimisations: Energy minimisation was carried out by 
using the AMBER 9 suite of programs in conjunction with the 
AMBER 99 forcefield. [28] The non-phosphorylated peptide, DWSFYL- 
LYYTEFT, containing an acetyl N terminus and a NH 2 C terminus, 
was first built as a monomeric strand; the steric zippers corre- 
sponding to class 7 and 8 were then constructed as two sheets of 
eight (3-strands each, separated by 12 and 4.8 A, respectively, by 
using the NAB molecular building tool to generate two antiparallel, 
shifted fibril models. [29] Side-chain C a dihedral angle values %^ were 
altered by using the Richardson rotamer library. [30] The minimisa- 
tions were carried out in implicit solvent with the generalised 
Born/Surface Area (GB/SA) method, [31] by using the Tsui and Case 
parameters/ 321 and an interaction cut-off of 25 A. The energy mini- 
misation used steepest descent for ten cycles, followed by a switch 
to conjugate gradient for the remaining 9990 cycles. Molecular 
representations showing structures at the end of the minimisation 
of the structures were produced by using PyMOL. [33] 

MD simulations: Trajectories were first obtained by using the GB/ 
SA implicit solvent model, as the absence of solvent friction from 
discrete water molecules will accelerate the in silico collapse of an 
ordered array if repulsive interactions outweigh favourable inter- 
actions between the peptide strands. 

The class 7 models obtained from the GB/SA minimisations were 
pre-equilibrated for 80 ps with coordinate, intersheet and hydro- 
gen bond restraints, then 1 ns without the coordinate restraints, 
and finally between 10 and 20 ns trajectories were obtained by 
running the structures unrestrained. A PDB file of the last snapshot 
was made and fully protonated phosphotyrosine residues (Y0P) 
were introduced at positions 5, 8 and 9 by using Amber parame- 
ters. [34] These new structures were minimised, pre-equilibrated with 
restraints and run for 10 ns in GB/SA unrestrained. A PDB file of 
each structure obtained was made. The PDB files of the Y0P fibrils 
were solvated by using a TIP3P solvation box, [35] minimised, pre- 
equilibrated and run unrestrained for 10 ns. 
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Fibrillar structures containing a mono- or di-anionic phosphate 
group (YIP and Y2P, respectively) were obtained by replacing Y0P 
in the PDB files from the GB/SA trajectories with either YIP or 
Y2P. [34] The structures generated were solvated with a TIP3P solva- 
tion box, [35] neutralised with counter-ions (Na + ), minimised in two 
steps, pre-equilibrated with coordinate restraints and finally run 
unrestrained for 10 ns. 



A cknowledgemen ts 

This work was funded by the Wellcome Trust (080706/Z/06/Z). 
The authors thank Dr. John White for his contribution towards 
the synthesis work, and Dr. Joshua Berryman for his help with 
the simulation work on fibril models. 

Keywords: molecular dynamics • peptides • phosphopeptides • 
protein folding • self-assembly 

[1] M. B. Pepys, Annu. Rev. Med. 2006, 57, 223-241. 
[2] M. Stridsberg, E. Wilander, Acta Oncologica 1991, 30, 451 -456. 
[3] G. G. Glenner, C. W. Wong, Biochem. Biophys. Res. Commun. 1984, 720, 
885-890. 

[4] K. A. Conway, J. D. Harper, P.T. Lansbury, Biochemistry 2000, 39, 2552- 
2563. 

[5] a) L C. Serpell, Biochim. Biophys. Acta 2000, 7502, 16-30; b) O. S. Makin, 

P. Sikorski, L. C. Serpell, Curr. Opin. Chem. Biol. 2006, 70, 417-422. 
[6] a) C. P. Jaroniec, C. E. MacPhee, N. S. Astrof, C. M. Dobson, R. G. Griffin, 

Proc. Natl. Acad. Sci. USA 2002, 99, 16748-16753; b) A. T. Petkova, Y. 

Ishii, J. J. Balbach, O. N. Antzutkin, R. D. Leapman, F. Delaglio, R. Tycko, 

Proc. Natl. Acad. Sci. USA 2002, 99, 16742-16747; c) R. Tycko, Curr. Opin. 

Chem. Biol. 2000, 4, 500-506. 
[7] J. L. Jimenez, J. L. Guijarro, E. Orlova, J. Zurdo, C. M. Dobson, M. Sunde, 

H. R. Saibil, EM BO J. 1999, 18, 815-821. 
[8] M. Torok, S. Milton, R. Kayed, P. Wu, T. Mclntire, C. G. Glabe, R. Langen, J. 

Biol. Chem. 2002, 277, 40810-40815. 
[9] a) R. Diaz-Avalos, C. Long, E. Fontano, M. Balbirnie, R. Grothe, D. Eisen- 

berg, D. L D. Caspar, J. Mol. Biol. 2003, 330, 1165-1175; b)0. S. Makin, 

E. Atkins, P. Sikorski, J. Johansson, L C. Serpell, Proc. Natl. Acad. Sci. USA 

2005, 702,315-320. 
[10] X.J. Lu, P. L. Wintrode, W. K. Surewicz, Proc. Natl. Acad. Sci. USA 2007, 

704, 1510-1515. 

[11] A. D. Williams, E. Portelius, I. Kheterpal, J.T. Guo, K. D. Cook, Y. Xu, R. 

Wetzel, J. Mol. Biol. 2004, 335, 833-842. 
[12] M. R. Sawaya, S. Sambashivan, R. Nelson, M. I. Ivanova, S. A. Sievers, M. I. 

Apostol, M.J. Thompson, M. Balbirnie, J. J. W. Wiltzius, H.T. McFarlane, 

A. O. Madsen, C. Riekel, D. Eisenberg, Nature 2007, 447, 453-457. 
[13] A. Schneider, J. Biernat, M. von Bergen, E. Mandelkow, E. M. Mandelkow, 

Biochemistry 1999, 38, 3549-3558. 
[14] K. E. Paleologou, A. W. Schmid, C. C. Rospigliosi, H. Y. Kim, G. R. Lamber- 

to, R. A. Fredenburg, P. T. Lansbury, C. O. Fernandez, D. Eliezer, M. 

Zweckstetter, H. A. Lashuel, J. Biol. Chem. 2008, 283, 16895-16905. 
[15] a) H. Kuhnle, H. G. Borner, Angew. Chem. 2009, 727, 6552-6556; Angew. 

Chem. Int. Ed. 2009, 48, 6431-6434; b) M. Broncel, J. A. Falenski, S. C. 

Wagner, C. P. R. Hackenberger, B. Koksch, Chem. Eur. J. 2010, 76, 7881 - 

7888; c) M. Broncel, S. C. Wagner, C. P. R. Hackenberger, B. Koksch, 

Chem. Commun. 2010, 46, 3080-3082. 
[16] S. Jones, J. Manning, N. M. Kad, S. E. Radford, J. Mol. Biol. 2003, 325, 

249-257. 

[17] a) A.M. Fernandez-Escamilla, F. Rousseau, J. Schymkowitz, L. Serrano, 

Nat. Biotechnol. 2004, 22, 1302-1306; b) F. Rousseau, J. Schymkowitz, L 

Serrano, Curr. Opin. Struct. Biol. 2006, 76, 118-126. 
[18] K. E. Routledge, G. G. Tartaglia, G. W. Piatt, M. Vendrusco, S. E. Radford, J. 

Mol. Biol. 2009, 389, 776-786. 
[19] a) L A. Carpino, G. Y. Han, J. Am. Chem. Soc. 1970, 92, 5748 - 5753; b) W. 

Chan, C. W. Chan, P. D. White, Fmoc Solid Phase Peptide Synthesis: A 

Practical Approach, Oxford University Press, 2000. 



280 www.chembiochem.org © 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim ChemBioChem 2012, 13, 271 -281 



Phosphorylation as a Probe of Fibril Architecture 



FULL PAPERS 



[20] S. B. H. Kent, Annu. Rev. Biochem. 1988, 57, 957-989. 

[21] a)C Hyde, T. Johnson, D. Owen, M. Quibell, R. C. Sheppard, Int. J. Pept. 
Protein Res. 1994, 43, 431-440; b)T. Johnson, M. Quibell, R. C. Shep- 
pard, J Pept Sci 1995, 7, 11 -25; c)W. R. Sampson, H. Patsiouras, N.J. 
Ede, J. Pept. Sci. 1999, 5, 403-409; d)W. G. Zeng, P.O. Regamey, K. 
Rose, Y. Z. Wang, E. Bayer, J. Pept. Res. 1997, 49, 273-279. 

[22] a) A. Baumketner, M. G. Krone, J.-E. Shea, Proc. Natl. Acad. Sci. USA 2008, 
105, 6027-6032; b) R. Tycko, K. L. Sciarretta, J. P. R. O. Orgel, S. C. Mere- 
dith, Biochemistry 2009, 48, 6072-6084; c) A.T. Petkova, R. D. Leapman, 
Z. Guo, W.-M. Yau, M. P. Mattson, R. Tycko, Science 2005, 307, 262-265; 
d) M.-K. Cho, H.-Y. Kim, C. O. Fernandez, S. Becker, M. Zweckstetter, Pro- 
tein Sci. 2011,20, 387-395. 

[23] Which symmetry relation leads to an in-register or shifted structure de- 
pends upon whether the peptide has an even or odd number of amino 
acid residues. Eisenberg's steric zipper classes are defined by symmetry 
but are exemplified by peptides containing an even number of resi- 
dues. This means that in Eisenberg's examples shifted antiparallel sheet 
arrangements occur in class 5 and 6 architectures. In our case the pep- 
tide contains an odd number of amino acid residues and consequently 
the shifted arrangements within sheets occur in architectures with 
class 7 and 8 symmetries. 

[24] E. Gazit, FASEB J. 2002, 76, 77-83. 

[25] a) A. Janer, A. Werner, J. Takahashi-Fujigasaki, A. Daret, H. Fujigasaki, K. 
Takada, C. Duyckaerts, A. Brice, A. Dejean, A. Sittler, Hum. Mol. Genet. 
2010, 79, 181 -195; b) L. A. Omtvedt, D. Bailey, D. V. Renouf, M.J. 
Davies, N. A. Paramonov, S. Haavik, G. Husby, K. Sletten, E. F. Hounsell, 
Amyloid 2000, 7, 227-244; c) K. E. Paleologou, A. Oueslati, G. Shakked, 

C. C. Rospigliosi, H. Y. Kim, G. R. Lamberto, C. O. Fernandez, A. Schmid, F. 
Chegini, W. P. Gai, D. Chiappe, M. Moniatte, B. L. Schneider, P. Aebischer, 

D. Eliezer, M. Zweckstetter, E. Masliah, H. A. Lashuel, J. Neurosci. 2010, 
30, 3184-3198. 



[26] a)Y. Miller, B. Ma, R. Nussinov, Proc. Natl. Acad. Sci. USA 2010, 707, 
9490-9495; b)Y. Miller, B. Ma, R. Nussinov, Chem. Rev. 2010, 770, 4820- 
4838; c)T. Eichner, A. P. Kalverda, G. S. Thompson, S. W. Homans, S. E. 
Radford, Mol. Cell 2011, 47, 161 -172. 

[27] a) C. Wasmer, A. Lange, H. Van Melckebeke, A. B. Siemer, R. Riek, B. H. 
Meier, Science 2008, 379, 1523-1526; b) A. K. Paravastu, R. D. Leapman, 
W.-M. Yau, R. Tycko, Proc. Natl. Acad. Sci. USA 2008, 705, 18349-18354. 

[28] D. A. Case, T. A. Darden, T. E. Cheatham III, C. L. Simmerling, J. Wang, 
R. E. Duke, R. Luo, K. M. Merz, D. A. Pearlman, M. Crowley, R. C. Walker, 
W. Zhang, B. Wang, S. Hayik, A. Roitberg, G. Seabra, K. F. Wong, F. Paesa- 
ni, X. Wu, S. R. Brozell, V. Tsui, H. Gohlke, L. Yang, C. Tan, J. Mongan, V. 
Hornak, G. Cui, P. Beroza, D. H. Mathews, C. Schafmeister, W. S. Ross, 
P. A. Kollman, AMBER, University of California, San Francisco, 2006. 

[29] Modeling Unusual Nucleic Acid Structures, T. J. Macke, D. A. Case, ACS 
Symp. Ser. 1998, 682, pp 379-393. 

[30] a) S. C. Lovell, I. W. Davis, W. B. Adrendall, P. I. W. de Bakker, J. M. Word, 
M. G. Prisant, J. S. Richardson, D. C. Richardson, Proteins Struct. Funct. 
Bioinf. 2003, 50, 437-450; b) S. C. Lovell, J. M. Word, J. S. Richardson, 
D. C. Richardson, Proteins Struct. Funct. Bioinf. 2000, 40, 389-408. 

[31] G. D. Hawkins, C.J. Cramer, D. G. Truhlar, J. Phys. Chem. 1996, 700, 
19824-19839. 

[32] V. Tsui, D. A. Case, Biopolymers 2000, 56, 275-291. 

[33] The PyMOL Molecular Graphics System, Version 1.3, Schrodinger LLC. 

[34] a) J. W. Craft, G. B. Legge, J. Biomol. NMR 2005, 33, 15-24; b) N. Homey- 
er, A. H.C. Horn, H. Lanig, H. Sticht, J. Mol. Model. 2006, 12, 281 -289. 

[35] W. L. Jorgensen, J. Chandrasekhar, J. D. Madura, R. W. Impey, M. L. Klein, 
J. Chem. Phys. 1983, 79, 926-935. 



Received: September 26, 2011 
Published online on December 15, 2011 



ChemBioChem 2012, 13, 271 -281 © 2012 Wiley-VCH Verlag GmbH&Co. KGaA, Weinheim www.chembiochem.org 281 



